May 27, 2009 13:15 WSPC/INSTRUCTION FILE schee-stuchlik 



Optical phenomena in the field of braneworld Kerr black holes 



Jan Schee 

Institute of Physics, Faculty of Philosophy and Science, Silesian University in Opava, 
Bezrucovo nam. 13, CZ-1^6 01 Opava, Czech Republic 

schee@email.cz 
Zdenek Stuchh'k 

Institute of Physics, Faculty of Philosophy and Science, Silesian University in Opava, 
Bezrucovo nam. 13, CZ-IJ^.6 01 Opava, Czech Republic 

Zdenek. Stuchlik@fpf. slu. cz 



We study the influence of the tidal charge parameter of the braneworld models onto 
some optical phenomena in rotating black hole spacetimes. The escape photon cones are 
determined for special families of locally non-rotating, circular geodetical and radially 
free falling observers. The silhuette of a rotating black hole, the shape of an equatorial 
thin accretion disk and time delay effect for direct and indirect images of a radiation hot 
spot orbiting the black hole are given and classified in terms of the black hole rotational 
and tidal parameters. It is shown that rising of negatively- valued tidal parameter, with 
rotational parameter fixed, generally strenghtens the relativistic effects and suppresses 
the rotation induced asymmetries in the optical phenomena. 
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1. Introduction 

One way of realizing theories describing gravity as a truly higher-dimensional in- 
teraction becoming effectively 4D at low-enough energies is represented by the 
braneworld models, where the observable universe is a 3-brane (domain wall) to 
which the standard model (non-gravitational) matter fields are confined, while grav- 
ity field enters the extra spatial dimensions the size of which may be much larger 
than the Planck length scale lp ~ 10 -33 cm El. 

As shown by Randall and Sundrum^H gravity can be localized near the brane at 
low energies even with a non-compact, infinite size extra dimension with the warped 
spacetime satisfying the 5D Einstein equations with negative cosmological constant. 
Then an arbitrary energy- momentum tensor could be allowed on the brane b ^ J . 

The Randall- Sundrum model gives 4D Einstein gravity in low energy limit, and 
the conventional potential of Newtonian gravity appears on the 3-brane with high 
accuracy SB. Significant deviations from the Einstein gravity occur at ver y high en- 
ergies, e.g., in the very early universe, and in vicinity of compact objects BEEnEC ] 
Gravitational collapse of matter trapped on the brane results in black holes mainly 
localized on the brane, but their horizon could be extended into the extra dimen- 
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sion. The high-energy effects produced by the gravitational collapse are disconnected 
from the outside space by the horizon, but they could have a signature on the brane, 
influencing properties of black holes ^4 There are high-energy effects of local char- 
acter influencing pressure in collapsing matter, and also non-local corrections of 
"backreaction" character arising from the influence of the Weyl curvature of the 
bulk space on the brane - the matter on the brane induces Weyl curvature in the 
bulk which makes influence on the structures on the brane due to the bulk gravi- 
ton stresses ^SJ The combination of high-energy (local) and bulk stress (non-local) 
effects alters significantly the matching problem on the brane, compared to the 4D 
Einstein gravity; f or sp herical objects, matching no longer leads to a Schwarzschild 
exterior in general I 2BHH1 The Weyl stresses induced by bulk gravitons imply that 
the matching conditions do not have unique solution on the brane; in fact, knowl- 
edge of the 5D Weyl tensor is needed as a minimum condition for uniqueness Q21@ 

1191 

Some solutions for spherically symmetric black holes and uniform density stars 
021 have been discussed. It has been shown that in the black hole case the matching 
conditions could be satisfied and the bulk effects on the black hole spacetimes could 
be represented by a single "brany" parameter. 

Assuming spherically symmetric metric induced on the 3-brane, the constrained 
effective gravitational field equations on the brane could be solved, giving Reissner- 
Nordstrom static black hole solutions endowed with a braneworld parameter b hav- 
ing character of a "tidal" charge, instead of the standard electric charge parameter 

Q 2 EI. The tidal charge 

can be both positive and negative, however, there are some 
indications that negative tidal charge should properly represent the "backreaction" 
effects of the bulk space Weyl tensor on the brane ^21 

The stationary and axisymmetric solutions describing rotating black holes lo- 
calized in the Randall- Sundrum braneworld were derived in ^ having the metric 
tensor of the Kerr-Newman form with a tidal charge describing the 5D correction 
term generated by the 5D Weyl tensor stresses. The tidal charge has an "electric" 
character again and arises due to the 5D gravitational coupling between the brane 
and the bulk, reflected on the brane through the "electric" part of the bulk Weyl 
tensor - , in analogy with the spherically symmetric case 021. 

When both the tidal and electric charge are present the black hole spacetime 
structure is much more complex and additional off-diagonal metric components g r ^, 
g rt are relevant along with the standard g^t component, due to the combination of 
the local bulk effects and the rotational dragging. This distorts the event horizon 
which becomes a stack of non-uniformly rotating null circles having different radii 
at fixed while going from the equatorial plane to the poles ^. The uniformly 
rotating horizon is recovered for the rotation parameter a small enough where Kerr- 
Newman form of the metric tensor is allowed describing charged and slowly rotating 
black holes In the absence of rotation, the metric tensor reduces to the Reissner- 
Nordstrom form with correction term of non-local origin^. 



a At present, no exact 5D solution in the braneworld model is known. 
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Here we restrict our attention to the Kerr-Newman type of solutions describing 
the braneworld rotating black holes with no electric charge, since in astrophysically 
relevant situations the electric charge of the black hole must be exactly zero, or very 
small 1221. Then the results obtained in analysing the behaviour of test particles and 
photons or test fields around the Kerr-Newman black holes could be used assuming 
both positive and negative values of the braneworld tidal parameter b (used instead 
of charge parameter Q 2 ). 

The information on the properties of strong gravitational fields in vicinity of 
compact objects, namely of black holes, is encoded into optical phenomena of dif- 
ferent kind that enable us to make estimates of the black hole parameters, including 
its tidal charge, when predictions of the theoretical models are confronted with the 
observed data. From this point of view, the spectral profiles of accretion discs around 
the black holes in galactic binaries, e.g., in microqu asars, are most p romising 
along with profiled spectral lines in the X-ray flux | 25 | 5 | 42 | 21 | 27 | 50| i m p Q rtant in- 
formation could also be obtained from the quasiperiodic oscillations observed in the 
X-ray flux of some low- mass black hole binaries of Galactic origin , some expected 
intermediate black hole sources 22] Q r those observed in Galactic nuclei EH. In the 
case of our Galaxy centre black hole Sgr A*, we could be able to measure the optical 
phenomena in more detailed form as compared with the other sources, since it is 
the nearest supermassive black hole with mass estimated to be ~ 4 x 1O 6 M 121, en- 
abling to measure the " silhuette " of the black hole and other subtle GR phenomena 

EEE] 

In the present paper, we give an introductory study of the tidal charge influence 
on the optical phenomena near a rotating black hole. We focus our attention to 
some characteristic phenomena in close vicinity of the black-hole horizon, where 
the effects of the tidal charge could be in principle of the same order as those of 
the black hole mass and spin, contrary to the case of weak lensing effects. The light 
escape cones are given for families of astrophysically interesting sources, namely 
in locally non-rotating frames, and frames related to circular geodetical motion 
and radially free-falling sources in section 4 1221 The silhuette of the black hole is 
determined in section 5. Images of the accretion discs are determined in section 6 
using the transfer-function method. In Section 7, time delay of hot spot radiation 
is determined for direct and indirect images assuming circular geodetical motion 
in close vicinity of the black hole horizon. In Section 8 relevance of some effects 
is estimated for the Galaxy centre Sgr A* supermassive black hole. Concluding 
remarks are presented in Section 9. 



2. Gravitational field equations on the brane 

In the 5D warped space models of Randall and Sundrum, involving a non-compact 
extra dim ensio n, the gravitational field equations in the bulk can be expressed in 
the form 
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G AB = k 2 [-Ag AB + 5( x )(-*9ab +T ab )\, (1) 

where the fundamental 5D Planck mass Mp enters via k 2 = Stt/M^ A is the brane 
tension, and A is the negative bulk cosmological constant. Denoting x — % A as the 
fifth dimension coordinate , x — determines location of the brane in the bulk 
space, at the point of Z2 symmetry; qab = Qab — tlatlb is the induced metric on 
the brane, with ua being the unit vector normal to the brane. 

The effective gravitational field equations induced on the brane are determined 
by the bulk field equations ([1]), the Gauss - Codazzi equations and the generalised 
matching Israel conditions with ^-symmetry. They can be expressed as modified 

standard Einstein's equations containing additional terms reflecting bulk effects 

1001 

onto the brane 



G [iv — ^Qybv k + k S^ V Efjy , (2) 

where k 2 = 87r/Mj>, with Mp being the braneworld Planck mass. The relations of 
the energy scales and cosmological constants are given in the form 




4.7T 



A- 



4.7T 



X 2 



(3) 



Local bulk effects on the matter are determined by the "squared energy- momentum" 

(4) 



tensor S^ v , that reads 



Spv = —TT^ - -T^ a T ua + —g/jLu {3T af3 T a(3 - T 2 ) , 



12" 4 _/i ' 24' 

while the non-local bulk effects are given by the tensor E^ v representing the bulk 
Weyl tensor Cabcd projected onto the brane, whereas 



£ab = C A BCDn c n D '. (5) 

Symmetries of the Weyl tensor imply that £[ab] = &a A = anc ^ £abti b = 
0. Therefore, on the brane, x ~^ 0? there is Sab —> £\lv^a"^b '• ^ ne tensor 
reflects influence of the non-local gravitational effects in the bulk, including the 
tidal ( "Coulomb ") and transverse traceless (gravitational wave) imprints of the free 
gravitational field of the bulk. 

We restrict our attention to the vacuum (at both bulk and brane) solutions of 
the gravitational field equations on the brane. Assuming zero cosmological constant 
on the brane (A = 0) we arrive to the condition 



A 



4^A 2 
~3Mj>" 



(6) 
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In the absence of matter fields, there is T^ v = = S^ V1 i.e., we are not interested 
in the properties of the squared energy-momentum S^ v representing local effects 
of the bulk. In the vacuum case, the effective gravitational field equations on the 
brane reduce to the formal 

iV = -V> v = ° = v ( ? ) 

implying divergence constraint ^^Hl 



V%„ = (8) 

where V M denotes the covariant derivative on the brane. 

The equation (|5J) represents Bianchi identities on the brane, i.e., an integrability 
condition for the field equations = — £^JM For stationary and axisymmetric 
(or static, spherically symmetric) solutions Eqs. (0) and (|5J) form a closed system 
of equations on the brane. 

The 4D general relativity energy- momentum tensor T^ v (with = 0) can be 
formally identified to the bulk Weyl term on the brane due to the correspondence 



k^T^ u <-> £fiv (9) 

The general relativity conservation law V^T^ V = then corresponds to the con- 
straints equation on the brane (|BJ). This behaviour indicates that Einstein-Maxwell 
solutions in general relativity should correspond to branew orld v acuum solutions. 
This was indeed shown in the case of Schwarzchild (R-N) I2SU2I and Kerr (K-N) 
spacetimes - . In both of these solutions the influence of the non-local gravitational 
effects of the bulk on the brane are represented by a single " braneworld" parameter 
b. The Coulomb-like behaviour in the Newtonian potential 



inspired the name tidal charge 



3. Null geodesies in Kerr spacetime with a tidal charge 
3.1. Geometry 

Following the work of and using the standard Boyer-Linquist coordinates (£, r, #, 
(/?), we can write the line element of Kerr black- hole (or naked singularity) spacetime 
on the three-brane in the form 

-2 2Mr ~ b sm 2 6dtd4>, (11) 
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where 



E = r 2 + a 2 cos 2 6 (12) 
A = r 2 - 2Mr + a 2 + b (13) 
A= (r 2 +a 2 ) 2 -a 2 Asin 2 #. (14) 

M is the mass parameter, a = J/M is the specific angular momentum and the 
braneworld prarameter b is the charge representing imprint of non-local grav- 
itational effects from the bulk space. The metric ([IT]) has the same form as the 
Kerr- Newman metric, where the tidal charge is replaced by the squared electric 
charge, Q 2 . The stress tensor on the brane E^ v takes the form^ 



E t t = -E^ = -^[E-2(r 2 + a 2 )}, (15) 

e; = —E g e = - A (i 6 ) 

V = -(r 2 + a 2 )sin 2 <9, (17) 
£/ = -Ay + a 2 )sin 2 

that is fully analogical (b — > Q 2 ) to the components of energy- momentum tensor 
for Kerr-Newman spacetimes in Einstein's general relativity^. 

The roots of A = identify the type of braneworld Kerr spacetime. There are 
two possibilities, a black hole or a naked singularity. By introducing a 2 /M 2 — > a 2 , 
b/M 2 — > 6, r + /M — > r_|_, or putting M = 1, we write the roots of A = in the form 



r + = 1 + \/l — <^ 2 — ^ 7 (outer horizon) (19) 

and 

r_ = 1 — yl — a 2 — 6, (inner horizon). (20) 

The metric given by the line element (jTTJ) determines the geometry of rotating black 
hole in braneworld universe if 



l>a 2 +6. (21) 

The strong inequality refers to the case of two horizonts r + and r_. For extreme 
black holes (1 = a 2 + b) the horizons coincide r+ = r_ = 1. 

It is clear that for 6 > the loci of the inner horizon r_ are always positive. But 
for b < 0, the loci of the inner horizon can also be at negative r, as illustrated in 
Figure [U 

Notice that a 2 > 1 is not allowed for standard black holes and for b > ^Sl, 
but such a possibility appears for b < 0. The rotational parameter of extreme black 
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Fig. 1. The plot of the inner horizont radius r_ as a function of tidal charge parameter b for 
three representative values of rotational parameter a 2 = 0.5, a 2 = 1.0 and a 2 = 1.5. 



holes is given by a 2 = 1 - b. The case of 1 < a 2 + b refers to the braneworld Kerr 
naked singularities. 

In this paper we focus on astrophysicaly interesting case of black holes, with 
emitting sources and observers located above the outer horizont. 



3.2. Carter's equations 

In order to study the optical effects in braneworld Kerr spacetimes, we have to solve 
equations of motion of photons given by the null geodesies of the spacetime under 
consideration. The geodesic equation reads 



where k^ = is the wave vector tangent to the null geodesic and w is the affine 
parameter. The normalization condition reads g^k^k" = 0. Since the components 
of the metric tensor do not depend on ip and t coordinates, the conjugate momenta 



K = g^ v k v = *, (23) 
h = gtuk u ee —E, (24) 

are the integrals of motion. Carter found another integral of motion K as a separa- 
tion constant when solving Hamilton- Jacobi equation 



(25) 
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where he assumed the action S in separated form 

S = -Et + + S r (r) + S e (0). (26) 
The equations of motion can be integrated and written separatelly in the form 

E t = ± " /W) > (27) 

AO 



E d^ = ± ^ W ^>> ( 28 ) 
^ = -5. + ^, (29) 



dw sin 2 A 
diu 



dt (r 2 +a 2 )P^ 

= -aP^ + t , (30) 



where 



R(r) =P 2 R - AK, (31) 

w ^"-'-(M- <32) 

P R (r)=E(r 2 + a 2 )-a<S>, (33) 
= a£sin 2 6>-<!>. (34) 
It is usefull to introduce integral of motion Q by the formula 

Q = K-(E-a(t)) 2 . (35) 

Its relevance comes from the fact that in the case of astrophysically most important 
motion in the equatorial plane (O = tt/2) there is Q = 0. 

3.3. Radial and latitudinal motion 

The photon motion (with fixed constants of motion E, <I>, Q) is allowed in regions 
where R(r;E,<$>,Q) > and W(9]E,$,Q) > 0. The conditions R(r;E,$,Q) = 
and W(6; E, <£, Q) = determine turning points of the radial and latitudinal mo- 
tion, respectively, giving boundaries of the region allowed for the motion. Detailed 
analysis of the ^-motion can be found in ^ while the radial motion was anal- 
ysed (with restrictions implied by the ^-motion) in SU and Here we extend this 
analysis to the case of b < 0. 

The radial and latitudinal Carter equations read 

^ 2 / dr \ _ 2 2 ai2 a r p n \ i 2i 



'o5/J = r -«A] 2 -A[£-2aA + a 2 ], (36) 

S 2 / d9 y = c + a 2 COS 2 g _ _aj_ (37) 

VdW sin 2 6> V ; 
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where we have introduced impact parameters 

A = ^, (38) 

£=i = ^^=4 + A 2 , (39) 

and rescaled the affine parameter by w' = Ew. We assume a > 0. 

The reality conditions (dr/di(/) 2 > and (d0/dw f ) 2 > lead to the restrictions 
on the impact parameter C 



Cmin — — ^max-) (^^) 

where 



(a\-2r + b) 2 9 rt 
>C ma . = ^ t + r 2 + 2r - 6, (41) 



and 



c . = i A 2 for |A| > a, 

^""1 2a|A| -a 2 for |A| < a. 1 ' 



The upper (lower) constraint, Cmax(£>min)i comes from the radial- motion 
(latitudinal-motion) reality condition. The properties of the photon motion are de- 
termined by the behaviour of the surface C max (r; A, a, 6), as given by (|4T|) . The 
extrema of the surface C max (giving spherical photon orbits) are determined by 

X = X+ = r -l±^, (43) 
a 

a(r - 1) y ' 

The values of C max at these extreme points are given by 



C max ( A+ ) = C+ = 2r 2 + o 2 , (45) 
2r(r 3 - 3r + 46) + a 2 (r + l) 2 

The character of the extrema follows from the sign of d 2 Cmax / dr 2 . One finds that 



<9 2 r 8r 2 



= — , for A = A+, (47) 



<9r 2 A 

, for A = A_. (48) 



d 2 C max 8r 2 8r 



9r 2 A (r - l) 2 
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Fig. 2. Left: classification of Kerr spacetime in braneworld universe according to the values of 
a 2 + 6, b and n ex t (the number of local extrema of the curves A±, which is also the number of 
circular photon orbits in the equatorial plane). The classification regions are: I) for a 2 + b < 1 
and n ex t = 2, II) for a 2 + b < 1 and n ex t = 4, III) a 2 + b > 1 and b < 1 and n ex t = 2, IV) for 
a 2 + b > 1 and b > 1 and n ex t = 2, V) for a 2 + b > 1 and n ex t = 0, VI) for a 2 + b > 1 and 6 < 1 



and n e 



4, VII) for a + b > 1 and 6 > 1 and n ex t = 4. Right: zoom of the area in the dashed 



rectangle of the left plot, to cover regions VI and VII. 



Clearly, there are only minima of Crnax along for A = A+, corresponding to unstable 
spherical orbits. 

Further, we have to determine where the restrictions given by the latitudinal 
motion C m i n meet the restrictions on the radial motion jC m ax- We find that C m ax = 
A 2 (for |A| > a) is fullfilled where 



A = A 4 



a(b-2r±r 2 VA) 
r 2 — 2r + b 



while 



2a|A| - a 2 (for |A| < a) is fullfilled where 



(49) 



A = A = — [4r - r 2 - 2b - a 2 + 2a/A(6- 2r)]. 



(50) 



The extreme points of curves X± , which are also the intersection points of these 
curves with A_, are determined by the equation 



/(r;a,6) r 4 



6r d 



(9 + 4b)r 2 - 4(36 + a 2 )r + 4b(b + a 2 ) = 0. 



(51) 



The equation /(r; a, 6) = determines loci of the photon equatorial circular orbits; 
in an implicit form the radii are given by the condition 



2 r 2 (r-3) 2 + 46(r 2 - 3r + 6) 



4(r-6) ' (52) 

The maxima of the curve A, which also determine the intersections of curves A and 
A_ are located on r satisfying the equation 
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2r 3 - (3 + b)r 2 + 2br + a 2 = 0. (53) 

The braneworld Kerr spacetimes can be classified due to the properties of the 
photon motion as determined by the behaviour of the functions A±, A±, A. The clas- 
sification is governed by their divergences (i.e., by existence of the horizons) and the 
number of local extrema determining equatorial photon circular orbits n ex . There 
exist seven classes of the braneworld Kerr spacetimes, with the criteria of separation 
being a 2 + 6 >1, 6 >1 and n ex . The classification is represented in Figure El There 
are two different classes of the black- hole spacetimes, differing by the presence of 
the photon circular orbits under the inner horizon. However, in the astrophysically 
relevant region outside the outer horizon, both the classes are of the same character, 
having two unstable equatorial photon circular orbits, one corotating (at r p hi) and 
the other counter-rotating (at r p h2 > fphi)- The tidal charge b introduces no qual- 
itatively new feature into the behaviour of photon motion in the Kerr spacetimes, 
but the quantitative impact of b < with high magnitude are quite relevant, as 
shown in next sections. All the braneworld Kerr black holes with tidal charge b < 
belong to the class II discussed in the case of standard Kerr-Newman spacetimes 
We illustrate in Figures 03|5] functions A± , A± and A for such a black hole spacetime 
with parameters a = 0.9 and b = —1.0. In this case typical for braneworld Kerr 
black hole with b < there exist ten significiant values of A as given in Figures [3]- 

® 

For each interval of A as determined by the sequence of A^ - Aj introduced in 
Figure [3l there exists a characteristic type of behaviour of the restricting "radial" 
function C m ax and its relation to the "latitudinal" restricting function C m i n . They 
can be found in and will not be repeated here. 

The allowed values of the impact parameter C lie between the limiting functions 
Cmi n and Cmax- If the minimum C™™ = £ maa; (r min , A ) of the limiting function 
Cmax is less than the value of the limiting function C m im an incoming photon 
(k r < 0) travelling from infinity will return back for all values of Co G [Cmin] Cmax}- 
If C™™x > £"min, the incoming photon (k r < 0) travelling from infinity returns back 
if its impact parameter Co satisfies the condition Co > C™^x an d is captured by the 
black hole if Co < C™™ x - The minimum C™™x determines (with the particular value 
of A) a photon spherical orbit, i.e., a sphere where photons move with r = const but 
with varying latitude (and, of course, varying tp). When the condition Co = Cmin is 
satisfied simultaneously, the spherical photon orbit is transformed to an equatorial 
photon circular orbit. Photons with Co = C™™x coming from distant regions or 
regions close to the black hole horizon will wind up around the photon sphere. 
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Fig. 3. The graphs of the X±, X± and A functions are plotted for representative values of the 
parameters a = 0.9 and b = —1.0. The two dashed rectangle areas labeled with numbers 1 and 2 
are zoomed in the following figures. The horizontal gray dashed lines represent special values of 
the impact parameter A, denoted according to the text as A^-.-Aj. 
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Fig. 4. Left figure is the zoom of dashed area labelled 1 in previous figure. Right figure is the 
zoom of dashed area labelled 2 in previous figure. The dashed rectangle area here is zoomed in the 
next figure. 
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Fig. 5. The zoom of the dashed rectangle area in previous figure. 
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4. Light escape cones 

The optical phenomena related to accretion processes in the field of rotating black 
holes could be efficiently studied by using the notion of light escape cones of local 
observers (sources) that determine which portion of radiation emitted by a source 
could escape to infinity and, complementary, which portion is trapped by the black 
hole 36 . Here we focus our attention to four families of observers (sources) that are 
of direct physical relevance. 



4.1. Local frames of stationary and free-falling observers 

We consider three families of stationary frames, namely LNRF (Locally Nonro- 
tatig Frame), SF (Static Frame) and GF± (Circular Geodesic Frame) and one non- 
stationary frame, namely RFF (Radially Falling Frame). The LNRF are of high- 
est physical importance since the physical phenomena take the simplest form when 
expressed in suc h fram es, because the rotational spacetime effects are maximally 
suppressed there EES The GF± are directly related to Keplerian accretion discs in 
the equatorial plane of the spacetime, both corotating and counterrotating, while 
RFF are related to free-falling spherical accretion. The SF are fixed relative to 
distant observers. The GF± and RFF are geodetical frames, while SF and LNRF 
are generally accelerated frames. 

The radial and latitudinal 1-forms of the three stationary frame tetrads are 
common for all three stationary cases and read 



^^ = {0,V / S7A,0,0}, (54) 

u w = {0,0,^,0}. (55) 

LNRF correspond to observers with $ = (zero angular momentum observers). 
Their time and azimuthal 1-forms read 



,W = 



A V 

0,0,0 , (56) 



A 

= J -SIlnrf\I ^ sin 0, 0, 0, \/^sin0V. (57) 



where 



_ a(2Mr-b) 
* 'lnrf — ^ yoo) 

is the angular velocity of LNRF as seen by observers at infinity. 

The tetrad of SF corresponding to observers with Q = ,i.e. static relative to 
observers at infinity, is given by the formulae 
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, 2r-6 „ a(2r - b) sin 2 I 

1 =— ,0,0,— f ' =}, (59) 

S ^/S 2 - (2r - 6)S) J 



^ = |o,0,0,^ s _ff_ 6) B mg|. (60) 

The GF± observers move along (/^-direction in the equatorial plane with velocity 
VGF±(+---corotating, -...counterrotating) relative to the LNRF and with angular 
velocity ft relative to the static observers at infinity given by [SK] 



n ± = ± 2 fzL- . (6i) 

r z ± ayr — b 



The velocity Vgf± is given by 

V GF± = ± (r 2 + f)YTa(2r- b )_ 

VA(r 2 ±aY) V ; 

where Y = \fr — b. The standard Lorentz transformation of the LNRF tetrad gives 
the tetrad of GF± in the form 

(t) f r*-2r + b±aY (r 2 + a 2 )Y =p a(2r — b) \ 
u± —< — ,U,U,=F 7^ J> (p6) 



where 



Z± = r^r 2 - 3r + 2b±2aY. (65) 

Note that the GF± family is restricted to the equatorial plane, while LNRF are 
defined at any 6. 

The RFF observers have velocity 

V RFF = {V^ r \ VM} (66) 



as measured in LNRF. The radially free-falling (or free-escaping) observers starting 
(finishing) at infinity move with = const. Using the results of ™l, we find the 
velocity components of the free-falling frames in the LNRF frames 



V {9) = 0, (68) 
V iv) = 0. (69) 
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Clearly, the free-falling (free-escaping) observers move only radially in the LNRF, in 
analogy to particles radially moving in static frames of the Schwarzchild spacetimes. 
For the radially free-falling sources, the tetrad components uo^ and uo^ coincide 
with those of the LNRF tetrad, while and are transformed. The local 
Lorentz transformation of the LNRF to the RFF± tetrad yields 



(70) 



w ± ( T7 V^va 7 ' ' /' (71) 

4 e) = {0,0,Vs,0}, (72) 
Jrt = J -n LNRF i M sin 0, 0, 0, \ sin \ . (73) 




(r) 1 /AS T ^ I Yj 



4.2. Construction of escape cones 




Fig. 6. Definition of directional angles ao, (3q and 70 in a local frame. Vectors e r , eg, e^ are the 
basic tetrad vectors. Position of the observer (source) is given by the coordinates (7*0, 0o)- Vector k 
represents a photon as observed by the observer in the given tetrad and vector k! is its projection 
into the plane ( eg, e^). 



For each direction of emission in the local frame of a source, there is a corre- 
sponding pair of values of the impact parameters A and C which can be related to 
the directional cosines of the photon trajectory in the local frame at the position 
of the source. Of course, the analysis of the turning points of the radial motion of 
photons, presented in the previous section, is crucial in determining the local escape 
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cones as the boundary of the escape cone is given by directional angles related to 
spherical photon orbits. 

Projection of a photon 4-momentum k onto the local tetrad of an observer is 
given by the formulae 

fcW = -k (t) = 1, (74) 

&( r ) = k( r ^ = cosao, (75) 

— fc(0) = sinaocos/3o, (76) 

k M = fc(y,) = sina sin/?o, (77) 

where ao, (3o are directional angles of the photon in the local frame (see Figure [6]) 
and COS70 = sinaosin/^o- In terms of the local tetrad components of the photon 
4-momentum and the related directional angles, the conserved quantities, namely, 
the azimutal momentum <1>, energy E and K read 



* = *V = -LjVykW + u {r \kW + + u^\k^\ (78) 

E = —k t = J*\kW - J r \kW - J e \kW - J t] ^\ (79) 

K = X ^ E ( r2 + a2) " a ^ ~ ^ k ^} ' ( 8 °) 
The impact parameters A and C defined by relations ([38]) and (|39|) are thus fully 
determined by any double, D, of angles from the set M = [qq? Ab 7o]- 

Having defined the source frame, we can construct light escape cones assuming 
fixed coordinates of the source r*o, #o- Their construction proceedes in the following 
steps: 

• for given D, say D = [ao, A)], we calculate A = A(ao, fa), 

• A determines the behaviour of C max = C max (r; A), 

• from the analysis presented in the previous section we calculate minimum 
of Cmax, which reads C™ a n x = C max (r min ; A), 

• we search for such a double D which satisfies equation Co(ao,/3o) = 

Cmaxymin'i A). 

Here, we present in detail the construction of light escape cones in particular case of 
the LNRF. The procedure is analogous for the other stationary frames and simply 
modified for the free-falling frames, being radius dependent. 

The impact parameter A expressed in terms of the angle 70, related to the 
LNRF, reads 

Ao = 1 , (81) 

l&LNRFO + Ao sin 0o cos 7 o 

where index '0' refers to the frame with coordinates [7*0, #o]- The minimum of C m ax 
is located at 
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7o 



Fig. 7. Left. The functions Cmax and £ m i n = Aq are plotted together with representative con- 
stant functions C\ and £2 to demonstrate the construction of the photon escape cone. Right. The 
intersections of Cmax(lo) with A 2 (7o) give the interval of relevant values of 70 G [imin'ilmax]- 



i ' -x for A < A ^ 

where 



_ J V o A - a 2 for A > A G ~ , N 

~ I 1 _ |i + & +v . \ / \ „ _ „ 



fei = a 2 + 26 + aA - 3, (83) 

k 2 = |27(1 - a 2 - 6) + 2>/3^27(l- a 2 -6) 2 + /^J 7 . (84) 

The relevant values of C lie between Cmax and £ m i n determined by Eqs flHJ) and 
(|^2|) . The intersections of functions C max = £ ma:E (7o) and £ m i n (7o) give the relevant 
interval of angles 7 G [7mm, 7maJ (see Figure [7J). For each 7 from [7m m ,7 ma J we 
calculate minimal value of the photon impact parameter C for which the photon 
reaches the turning point r m i n and escapes to infinity. This minimal value is the 
minimum of C max which is located at r min , eg. C max = C max (r m i n ; A (7o), a, &), 
where r m i n is given by ([82]) . Now we can calculate the value of ao using equation 

cosa ° = I« = ^) 7Z- (85) 

We arrive to the formula 



onan , /T~ v/(rg+ffl 2 -fflAo)2-Ao(£^" x -2fflA + ffl^) 

cos c^o = ±y/io 9 9 9 9 o ' \*®) 

-a(a sin #0 — Ao)Ao + (rfi + a )(?o + a — «Ao^ 

where A = A(r ,# ), A = A(r ) and C™™ x = C max (rmin] Ao, a, b). The angle /? 
can be calculated from the formula ([77]) . In this way we obtain angles from the 
arc (3q G (— 7r/2; tt/2). The remaining arc (3q G (7r/2;37r/2) can be obtained by 
turning the arc /?o £ ( — 7r /2;7r/2) around the symmetry axis determined by angles 
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/?o = — 7r/2 and /3o = tt/2. This procedure can be done because photons released 
under angles (3q and 7r — (3q have the same constants of motion. Clearly, for sources 
under the radius corresponding to the corotating equatorial photon circular orbit, 
only outward directed photons with no turning point of the r-motion can escape. 
With radius of the source approaching the event horizon (r*o — ► r+), the escape 
cone shrinks to infinitesimal extension, except the case of extreme black hole 
For the other frames considered here, the procedure of the related light escape 
cone construction can be directly repeated, but with the relevant tetrad 1-form 
components being used in the procedure. 

In order to reflect properly the effect of the tidal charge b on the escape cone 
structure, we shall give the cones for black hole sequences of two kind: first we keep 
the spin a fixed and change 6, second we keep fixed "distance" to the extreme black 
hole states, i.e., a 2 + b is fixed, and both a and b are changed. The positive tidal 
charges have tendency to slightly increase the asymmetry of the cones as compared 
with 6 = case, keeping its character similar to the case of Kerr black holes (see 
next section). Therefore, we focus our attention to the influence of negative tidal 
charges. 





Fig. 8. Light escape cones as seen by LNRF in the vicinity of the braneworld kerr black hole. Top 
set of images is plotted for radial coordinate of emitter r e = 6M and bottom set for r e = 20 M. The 
rotational parameter a = 0.9981 is fixed and the representative values of the braneworld parameter 
b are (left), —1 (middle) and —3 (right). The shaded area represents photons captured by black 
hole. 



Behaviour of the LNRF escape cones in dependence on the braneworld parame- 
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Fig. 9. Light escape cones as seen by LNRF in the vicinity of the extreme braneworld kerr 
black hole. Top set of images is plotted for radial coordinate of emitter r e = 6M and bottom set 
for r e = 20M. The representative rotational and braneworld parameters [a 2 , b] are [1.0,0.0] (left), 
[2.0, — 1.0] (middle) and [4.0,— 3.0] (right). The shaded area represents photons captured by black 
hole. 

ter b (and the spin a) is represented in Figures [8] and [9l The complementary trapped 
cones, corresponding to photons captured by the black hole, are shaded. 

At a fixed radius expressed in units of M the extension of the trapped cone 
grows with descending of b to higher negative values and fixed spin a and mass M, 
demonstrating thus the growing gravitational pull of the black hole due to growing 
magnitude of the negative braneworld parameter. The same statement holds also 
in the case of extreme Kerr black holes, when a grows and b descends, while M 
is fixed. Clearly, the positive braneworld parameters have tendency to increase the 
asymmetry of the cones, while the negative ones symmetrize the escape cones with 
growing of On the other hand, the asymmetry of the escape cone grows with 
descending of b for extreme black holes (Figure [9]). 



Table 1. Table of relevant values of r ms and used 
in plots on Figs [10] and [TT] 



(a 2 , b) 


(0.9981,0.0) 


(0.9981,-1.0) 


(0.9981,-3.0) 


T'ms 


1.24M 


3.91M 


6.27M 




1.062M 


2.002M 


2.73M 
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Fig. 10. Escape cones of GF+ observers. Top images are plotted for observer (emitter) at r = r ms , 
middle images r = 10M and bottom images for r = 10 • r^. The value of a = 0.9981 is kept fixed. 
The representative values of b are (from left to right) 0.0, —1.0 and —3. 



Further, we represent the influence of the braneworld parameter on the escape 
cones for the circular (corotating) geodesic frames in Figure [101 Assuming astro- 
physically relevant sources in Keplerian accretion discs, their orbits must be lo cated 
above the marginally stable orbit r ms , determined implicitly by the condition 



4(r - 6) 3 / 2 =F r^/3r 2 - 2r(l + 2b) + 3b 
a = a ms (r- b) = _ . (87) 

Therefore, we construct the escape cones for observers at r = r ms (a, b) and at fixed 
radii. In the sequence of black holes with fixed spin a = 0.9981 (Figure fTO]) we include 
also a subsequence of escape cones constructed at the same relative distance from 
the black hole horizon in order to better illustrate the role of the tidal charge b. In 
the sequence of near-extreme black holes with a 2 + b = 0.9999 (Figure [TT]) the third 
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Fig. 11. Escape cones of GF+ observers. Top images are plotted for observer (emitter) at r = r m s 
and bottom images for r = 10M. The value of a 2 + b = 0.9999 is kept fixed. The representative 
values of (a 2 ; b) are (from left to right) (0.9999; 0.0), (1.9999; -1.0) and (3.9999; -3). 

sequence is not necessary as the black hole horizon is fixed at Th = 1.0 1M. Figures 
[TOl and [TT] demonstrate that the trapped cone expands as the tidal charge descends 
to lower negative values, both for black holes with fixed spin a and for near-extreme 
holes. On the other hand, considering the cones at r ms we can conclude that the 
descending tidal charge (b < 0) symmetrizes their shape for fixed a, but makes 
them strongly asymmetric for near-extreme black holes shrinking them strongly in 
the direction of the black hole rotation. 

Finally we demonstrate the relevance of the tidal charge b in the character 
of escape cones of the RFF- (comparing them with those related to LNRF) in 
Figure [T2j We construct the escape cones for two typical values of the tidal charge 
(b = 0, b = —3) in a sequence of radii where the free-falling source is radiating, 
demonstrating thus the combined growing influence of the black hole gravitational 
pull on the photon motion and the velocity of the free-falling source. In order to 
illustrate the phenomena in a clear way, we compare the RFF- escape cones to the 
corresponding LNRF escape cones. Clearly, the tidal charge descending to higher 
negative values makes stronger squeezing of the free-falling cones relative to the 
LNRF escape cones at any fixed radius. Notice that both the RFF- and LNRF 
cones are shifted due to the black hole rotational dragging. We again observe the 
tendency of negative brany parameters to symmetrize and squeeze the escape cones. 
At a fixed r, the escape cones become smaller for growing \b\ due to stronger gravity. 
For completeness we present sequence of both the RFF- and LNRF escape cones 
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at the three fixed radii for an extreme black hole with b — — 3 and a 2 = 4. We 
observe that both the RFF_ and LNRF cones are strongly shifted in the sense of 
the black hole rotation in vicinity of black hole horizon due to growing influence of 
the spin. The symmetrizing effect of descending values of negative tidal charge is 
canceled by strong influence of the rotational effects due to growing black hole spin. 
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Fig. 12. Comparison of the effect of the tidal charge b on the shape of light escape cones of locally 
nonrotating (dashed curves) frames and free falling (solid curves) frames. In the left column light 
escape cones are plotted for the tidal charge parameter 6 = and in the middle one the light 
escape cones are plotted for b = —3. The spin a = 0.9981 is kept fixed in both columns. The right 
column gives the sequence of the escape cones for an extreme black hole with [a 2 = 4; b = —3]. 
Emitting sources in all plots are moving in the equatorial plane. The radial distances of emitter 
are r e = 5M (top row), r e = 10M (middle row) and r e = 15M (bottom row). 



May 27, 2009 13:15 WSPC/INSTRUCTION FILE schee-stuchlik 



24 Jan Schee and Zdenek Stuchlik 

5. Silhuette of braneworld Kerr black hole 

In principle, it is of astrophysical importance to consider a black hole in front of 
a source of illumination whose angular size is large compared with the angular 
size of the black hole A distant observer will see a silhuette of the black hole, 
i.e., a black hole in the larger bright source. The rim of the black hole silhuette 
corresponds to photon trajectories spiralling around the black hole many times 
before they reach the observer. Of course, the shape of the silhuette enables, in 
principle, determination of the black hole parameters. But we have to be aware of 
the strong dependency of the silhuette shape on the observer viewing angle; clearly, 
the shape will be circular for observers on the black hole rotation axis, and its 
deformation grows with observer approaching the equatorial plane. 

Assuming that distant observers measure photon directions relative to the sym- 
metry center of the gravitational field, the component of the angular displacement 
perpendicular to the symmetry axis is given by —p^) /ptf) (for black hole rotating 
anticlockwise relative to distant observers), while for angular displacement parallel 
to the axis it is given by p^ /p^\ These angles are proportional to l/r*o, therefore, 
it is convenient to use the impact parameters in the form independent of ro ^ 

p&) A 

^ = -^0-77^ = — (88) 
pW sin #o 

and 

(3 = = \q + a 2 cos 2 - A 2 cot 2 ] 

p[t) L J 



C + a 2 cos 2 • 



A 2 



.1/2 



sin 2 #o 



(89) 



Photon trajectories reaching the observer are represented by points in the (a — (3) 
plane representing a small portion of the celestial sphere of the observer. 

The shape of the black hole silhuette is the boundary of the no-turning-point 
region, i.e., it is the curve C = C™™xW expressed in the (a — (3) plane of the 
impact parameters. For observers in the equatorial plane (0q = 7r/2), a = —A, 
P=(£-\ 2 ) 1 / 2 =q 1 / 2 . 

We consider that the black hole is observed by static distant observers. Therefore, 
we shall use the static frames introduced above. The silhuette of the black hole is 
quite naturally related to their trapped (escape) light cones. 

The marginal values of impact parameters Ao and £o(resp go) are obtained from 
the light escape cone. Using the stationarity of the braneworld Kerr spacetime we 
"shoot out" virtual photons from observer (static frame at very large distance ro) 
and we are looking for the light escape cone of this virtual source (using the results 
of the previous section). The trapped light cone of this virtual source is constructed 
from the light escape cone of the virtual source by transformations of directional 
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Fig. 13. Left figure. The (o:o,A)) plots of the silhuettes of braneworld Kerr black hole on a 
bright background for rotational parameter a? = 0.6 and four representative values of tidal charge 
parameter b = —3.0, b = —0.4, b = 0.0 and b = 0.4. The observer is located at ro = 10 4 M and 
#0 = 90°. Right figure. The silhuettes of extreme black holes for three representative values of 
braneworld parameter 6 = (solid), b = — 1 (dashed) and b = — 3 (dotted). Static observer is in 
equatorial plane at radial distance from the centre ro = 10 4 M. 

angle to ao = tt — ao and directional angle (3q to /3q = (3q. In this way we get 
marginal directions for received photons from bright background behind the black 
hole. Then we can use the formulas ([?8]) , ([79|) and ([80]) to calculate the marginal 
values of Ao and qo(£o) in order to obtain the silhuette of the braneworld Kerr 
black hole in the plane (a — /?), i.e., the set of doubles (ao,/3o) from equations 
([88]) and ([89]) . Here we plotted the silhuette directly from the trapped light cone 
(ao,/3o) on the observer's sky (ao sin/3o, <^o cos/?o). Note that the angle ao is the 
radial coordinate and the angle /3o is the polar coordinate in the polar graph of the 
silhuette. 

We shall give the silhuette of the black hole for observers located at fixed radius 
ro = 10 4 M that corresponds to the angular size of a ~ 1.4arcsec; for higher distances 
the angular size falls accordingly to the 1/ro dependence. 

First, we give an illustrative picture of the tidal charge influence on the silhuette 
properties for maximal inclination angle 0o = 90° when the black hole rotational 
effects are strongest (Figure [13]). We present a sequence of silhuettes for fixed black 
hole spin and varying b (left) and for extreme black holes with a 2 + b = 1 and 
both a, b varying (right). We clearly see that the positive tidal charge squeezes 
magnitude of the silhuette making its shape more asymmetric, while negative tidal 
charge enlarges silhuette's diameter symmetrizing its shape when a is fixed. For 
extreme black holes the silhuette asymmetry is kept but its extension grows with b 
descending to higher negative values. 

Second, there is a crucial effect of the viewing angle #o onto the shape of the 
black hole silhuette, demonstrated in Figure [14] for representative values of b and 
fixed spin a, and in Figure [15] for extreme black holes with parameters [a 2 = 1;6 = 0] 
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0° 0° 



Fig. 14. The silhuettes of rotating braneworld black hole on a bright background. Each image 
contains three black hole shapes for three representative values of observer's inclination angle 6q = 
{0° (solid), 45° (dashed), 90° (dotted)}, observer's radial coordinate tq = 10 4 M and the rotational 
parameter a 2 = 0.8. Top left image: b = 0.2. Top right image: b = 0.0. Bottom left image: b = —1.0. 
Bottom right image: b = —3.0. 

and [a 2 4;6 = -3]. 

The rotational effect on the shape of the silhuette grows with inclination angle 
growing and becomes strongest when 6q = tt/2; then the suppressing effect of the 
braneworld parameter is given in the most explicit form as demonstrated in Figure 

CE3 

The negative values of the braneworld parameter have the tendency to make 
the silhuette of a Kerr black hole (with a 2 fixed and for 7*0, 0q fixed) spherical, 
suppressing thus the rotational effects. However the symmetrizing effect of the tidal 
charge could be masked by symmetrizing effect of the viewing angle 6$. Therefore, 
it is very important for black hole parameter estimates to have observational limits 
on the value of 6q. 
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Fig. 15. The silhuettes of extreme rotating braneworld black holes on a bright background. Each 
image contains three black hole shapes for three representative values of observer's inclination angle 
#0 = {0° (solid), 45° (dashed), 90° (dotted)}, observer's radial coordinate ro = 10 4 M. Silhuettes on 
the left figure are plotted for extreme black holes with a 2 = 1 and 6 = and on the right side for 
a 2 = 4 and b = —3. 



180° 




magnitude of distorsion of Kerr black hole silhuette in braneworld universe. 



In order to characterize the influence of the tidal charge on the silhuette of a Kerr 
black hole we define two quantities in principle measurable by distant observers. The 
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shift s of the silhuette 

s a(Pm) sin(/3 m - tt), (90) 

and its ellipticity e 

_ ajp = 90°) +a((3 = 270°) 
' 2&(/3 m ) cos(/? m - tt) ' 1 j 

where /? m is defined by a(/3 m ) sm(/3 m — 7r) > a(/3) sin(/3 — 7r), V/? £ [7r/2,3/27r] 
i.e., it defines maximal extension of the silhuette in the x-direction. The definition 
of shift s and elipticity e is illustrated in Figure [T6l 

We calculated shift s and ellipticity e as functions of tidal parameter b for the 
Kerr black hole with rotational parameter a 2 = 0.9995 (see Figure [T7|). Clearly, 
these are quantities that could be measured and used for a black hole parameters 
estimates, if observational techniques could be developed to the level enabling the 
silhuette detailed measuring. We shall discuss such a possibility for the case of the 
supermassive black hole predicted in the Galaxy Centre (Sgr A*). 




Fig. 17. Top row. Left figure: the shift s = s(b) as a function of braneworld parameter b. Right 
figure: the ellipticity e = e(b) as a function of b. There are two curves on each image, one for 
observer inclination angle #o = 45° and second for #o = 89.9°. The rotational parameter of black 
hole is fixed to value a = 0.9995 and the radial coordinate of observer if ro = 10 4 M. Bottom 
row. The ellipticity e (left) and shift s (right) of the extreme black hole silhuette as functions of 
braneworld parameter b. Observer's coordinates are 0q = tt/2 and ro = 10 4 M. 
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6. Direct and indirect images of radiating disc 

Modelling of spectral line profiles of a thin radiating ring rotating in the equatorial 
plane of a braneworld Kerr black hole or light curve of an isotropically emitting 
point source orbiting such a black hole will give us information about the influence 
of the braneworld parameter b on the optical phenomena in the strong field regime 
1351. Here we restrict our attention to images of radiating discs. We can then, at 
least in principle, obtain estimates on the astrophysically acceptable values of the 
braneworld parameter b. 

6.1. Images of isoradial geodesies 

Calculating images of an accretion disc (ring) in the equatorial plane of a braneworld 
Kerr black hole is the first step to calculate the optical phenomena. Generally one 
could obtain a direct and an indirect image (see Figures fTHl and [T9|) , but in special 
cases the situation can be much more complicated due to complex character of the 
latitudinal and azimuthal photon motion. Here we focus our attention to the direct 
and indirect images of isoradial geodesies. 

In order to find all relevant positions of points forming the rotating ring on 
observer's sky, we have to find photon trajectories between the ring particles and 
the observer, i.e., we seek for such doubles of local observational angles [ao, Po] that 
satisfy the condition 
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Fig. 18. Direct image of the rotating ring in the equatorial plane at r e = 6M around braneworld 
Kerr black hole with rotational parameter a 2 = 0.5 for four representative values of tidal charge 
parameter b = —3.0, b = —0.4, b = 0.0 and b = 0.4. The observer is located at ro = 10 4 M and 
O = 85°. 



Iu(oLo,Po',nu,u a gn) ~ lM(oLo,0o',n,p,s) = 0. (92) 
Here we introduced the modified radial coordinate u = 1/r and cosine of latitudinal 

I O Q I 

coordinate fi = cos 0™. In the condition ([92|) n u is the number of turning points in u 
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Fig. 19. Indirect image of the rotating ring in the equatorial plane at r e = 6M around braneworld 
Kerr black hole with rotational parameter a 2 = 0.5 for four representative values of tidal charge 
parameter b = —3.0, b = —0.4, b = 0.0 and b = 0.4. The observer is located at ro = 10 4 M and 
O = 85°. 

coordinate, n is the number of turning points passed in \i coordinate, p = mod(n, 2), 
s = (1 — jjb sgn )/2. In terms of u and fi we define the functions I\j and Im by 



l[/(tt ,A;^,v) = } 1 &gn 



-u sgn + £;) for n u = 1 
u S gn C e for n u = 



(93) 



and 



Im (<*o , A) ; ^, s) = Ms^n 



/ +(-l) n+1 / + 



(94) 



+ (-l) s [(l - p)n +p[(l - s)(n - l) + «(n + 1 



J /x_ 



with 



du 



U(u) = 1 + (a 2 - A 2 - <z)u 2 + 2[(A 2 - a 2 ) 2 + <z]w 3 - 
- [q(a 2 +b) + b(a - \) 2 }u 4 



(95) 
(96) 



and 
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f - f -^=, (97) 
M(/i) = g + (a 2 - A 2 - ^)/i 2 - a 2 /A (98) 

6.2. Integration of photon trajectories 

We express the integrals ([95]) and ([97|) in the form of the standard elliptic integrals of 
the first kind. Rauch and Blandford presented the tables of reductions of ^-integrals 
and //-integrals for the case of photons in Kerr geometry ^1 Here we extended those 
reductions for the case of nonzero braneworld parameter b. Because the integration 
of the //-integral does not depend on braneworld parameter 6, the transformations 

IOOI 

are the same as in the case of Kerr metric 1 ^, but we include them for completeness. 

There are two cases we distinguish in latitudinal integral (see table [2]). In the 
first case there is one positive, M + > 0, and one negative, M_ < root of M(m 2 ) it 
implies that there are two turning points located symm etrically about the equatorial 
plane given by ±y/M+ (so called orbital motion I 8 | 13 | j n ^ e seGO nd case there are 
two positive roots, < M_ < M + of M(m 2 ), which implies that the latitudinal 
motion is constrained to the region above or below of the equatorial plane (so called 
vortical motion). The relevant reductions of the integral f™ dm'/ M{m') = 1m 
are stored in the table O 

For distant observers we distinguish five relevant cases of the radial integral. 
These cases depend on the character of roots of the quartic equation U{u) = 0, i.e., 
on the number of turning points {n u = or n u = 1) in the radial motion and the 
value of parameter q = q(a 2 +b) + b(a — A) 2 . We have arranged those transformations 
into table 03 

Denoting roots of the quartic equation U(u) = by /?2, /?3 and /?4, the 
meaning of each of the five cases is the following: 

• The case I: four distinct real roots of U(u) = forming the sequence 
Pi > fh > fh > and /?4 < 0. The value of modified constant of motion 
q > 0. 

• The case II: four real roots as in the case I but their values form the 
following order: (3\ > fa > and (3^ < /3s < 0. The value of modified 
constant of motion q < 0. 

• The case III: two real and two complex roots of U(u) = 0: (3\ being a 
complex root, fa — Pi and (3^ < fa < 0. The value of modified constant of 
motion q < 0. 

• The case IV: only complex roots: fa — Pi and P4 = 0s. The value of 
modified constant of motion q < 0. 

• The case V: two real and two complex roots of U(u) = 0: Pi > 0, P4 < 0, 
P2 being a complex root and Ps = 02- 
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Table 2. The reductions of dm' /^M(m') = I M 



Case tan ^ 



ci 



mi 



M_ < 



M_ > 



M + -M_ v /a 2 (M + -M_) 



M + — m 2 M_|_ — M_ i 
m 2 — M_ Mi a 2 " 



/M-L 



Table 3. The 



reductions of du' / y/U(u f ) = ly 



Case tan ^ 



ci 



II 



III 



IV 



V 



/(/3i 


-/3 3 )(u- 


04 ) 


(/3l 


-/3 4 )(/3 3 


-u) 






03 ) 



V (/3l-/3 3 )(/3 2 -n) 

2c 2 (u) 
|l-c 2 (u)| 

ii — c 3 

3(/3l)(l + c^) + c 2 (w-c 3 ) 
2c 2 (tQ 



\l-4(u)\ 



(/3i-/3 2 )(/3 3 -/3 4 ) 
(/3i-/3 3 )(/3 2 -/3 4 ) 

(/3 2 -/3 3 )(/3i-/3 4 ) 
(P1-P2XP4-P3) 

4c 4 c 5 -(/3 3 -/3 4 ) 2 -c 4 c 5 
4c 4 c 5 



/ c 4 -c 5 y 

V c 4 +c 5 / 



( C4+ c 5 ) 2 -(/3 1 -/3 4 )^ 
4c 4 c 5 



v /g(bl-63)(b2-64) 



v /-g(bl-b2)(b3-b4)] 



V-gc 4 c 5 



(c 4 + c 5 )^^g 



V<?c 4 c 5 



£4 
ft 
ft 
C3 
ft 



Table 4. Definitions for Table [3] 



Case 1 C2 



c 3 



III 



IV 



c5(u-(3 3 ) l 1/2 
c 4 (u-(3 4 ) J 



f 4[Q(/3 1 )] 2 -(c 4 -c 5 ) 2 ] 
I (c 4 +c 5 ) 2 -4[G(/3 1 )] 2 J 

f c4(^-/3 4 ) l 1 / 2 
Lc 5 (/3i-n)J 



2 >, 1/2 



^(ft) + C 2 3(ft) 



Using presented transformations we can write the integrals (|95|) and (|97j) in the 
form 



and 



r 1 



zdii = ci^ r (^; m) 



(99) 



(100) 



where T is the elliptic integral of the first kind and wi(resp depends on the 
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Table 5. Definitions for Tableland Table 0] 



Case 1 C4 



III {[5K(/3i)-/3 3 ] 2 + [S(/3i)] 2 } 1/2 {[»(/?!) -/34] 2 + [9G9i)] 2 } 1/2 

IV |[R(ft) - 5R(/3 3 )] 2 + [3(A) + 3(fe)] 2 } 1/2 {Wft) - 5R(fe)] 2 + [9(/3i) - 3(A)] 2 } 1/2 

V |[sR( / 3 2 )- /3l ]2 + [a( / 3 2 )]2} 1/2 {[5R(/3 2 )-/3 4 ] 2 + [9(/3 2 )] 2 } 1/2 
1 T/ie symbols dt(x) and 5s (x) refer to real and imaginary part of x here. 

case of root distribution of quartic equation U(u) = (resp. M(/i) = 0) as given in 
Table [3] (resp[2j). If, in the cases III and V, the value of 1 — c^u) < 0, we have to 
take instead of ([99]) the form 

f - 7 L=du = Cl (2/C(m) - ^(*; m)), (101) 

where K is the complete elliptic integral of the first kind. In the case that sign(/il • 
/i) < we have to take instead of ()100p the form 



/ 



/—=d/i = Cl (2/C(m) - .F(*; m)), (102) 

where m and ci are taken from table [21 We consider two basic possibilities of 
trajectories, namely those corresponding to direct and indirect images (Figures fT8l 
and [19]). 



6.3. Disc images 

It is very important to demostrate the influence of the braneworld parameter on 
the shape of images of rings in the equatorial plane representing parts of Keplerian 
accretion discs. Of course, as well known from the Kerr (and even Schwarzchild) 
black holes, the images strongly depend on the latitude of the observer. We calculate 
the direct and indirect images of flat discs and combined, full image of the disc for 
two representative values of viewing angle #o and appropriately chosen extension 
of radiating disc area. 

We include the effect of frequency shift into the calculated images of part of the 
Keplerian discs assumed to be radiating at a given fixed frequency. The frequency 
shift g is determined by the ratio of observed (Eq) to emitted (E e ) photon energy 

»=§=f4 < io3 > 

where Uq (u%) are components of the observer (emitter) 4- velocity and &o/x(^e/x) are 
components of the photon 4-momentum taken at the moment of emission (obser- 
vation). For distant observers Uq = (1,0,0,0). The emitter follows an equatorial 
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circular geodesies at r = r e , (9 e = tt/2. Therefore, u% = (i4, 0, 0, uf), with compo- 
nents given by 



ul 



1-1/2 

(104) 



_l(l_ a0) 2_ (r 2 +a 2 )0 2 + ^ (1 _, 

< = (105) 



where = dy/di is the Kepler ian angular velocity of the emitter related to distant 
observers, given by equation (|6ip . 

The frequency shift including all relativistic effects is then given by 

1 1/2 

1 - f (1 - atl) 2 - (r\ + a 2 )n 2 + £(1 - 2att) 

• — T^xn ' — ^ 

where A = —k^/k t is the impact parameter of the photon being a motion constant 
for an individual photon radiated at a specific position of the radiating disc; notice 
that g is independent of the second photon motion constant (impact parameter) 
q. Of course, depending on the position of the emitter along the circular orbit, 
the impact parameters A, q of photons reaching a fixed distant observer will vary 
periodically (see eg., ^). For each position of the emitter the impact parameters are 
determined by the procedure of integration of photon trajectories. 

The influence of the frequency shift in the disc images is demonstrated in Figures 
l20l and EH The role of the braneworld parameter is illustrated both for small (0q — 
30°) and high (0q = 80°) inclination angles. We consider two cases of the radiating 
disc extension: first one with fixed inner and outer radii, independent of the black 
hole parameters, and the second one when the inner radius is identified with the 
marginally stable orbits, depending on the black hole parameters. 

In order to map the frequency shift g into color palete we define modified fre- 
quency shift g = (g - g m in) / (gmax ~ gmin) where g min (g m ax) is the minimal (max- 
imal) value of frequency shift, which is fixed in a particular set of images. 

We can see from Figs. l20l and [2T1 that the negative tidal charge has the tendency 
to enlarge and symmetrize the disc images. 
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O 






Fig. 20. Radiating Keplerian disc images with fixed inner and outer radii. The modified frequency 
shift g = (g - g m in) I '{gmax - 9min), with gmin = 0.4 and gmax = 1-5, of the radiation emitted 
from the thin disk with inner radius r^ n = 7M and outer radius r ou t = 15M, encoded into colors is 
plotted for representative values of tidal charge parameter b = —3.0, 0.0 and inclination of observer 
#0 = 30°, 80°. In the left column direct images are ploted, the indirect images are ploted in the 
central column and the composition of direct and indirect images is plotted in the right column. 
The first two rows of images are plotted for the observer inclination 0q = 30° and the second two 
rows of images are plotted for the observer inclination 0o = 80°. Top row images are plotted for 
b = 0.0, the second row images are plotted for b = —3.0, the third row images are plotted for 
b = 0.0 and bottom row images are plotted for b = —3.0. 
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Fig. 21. Radiating Keplerian disc images with ri n = r ms . The modified frequency shift g = 
(g - 9min)/(9max - g m in), with g m i n = 0.2 and gmax = 1.8, of the radiation emitted from the 
thin disk with inner radius r^ n = r ms (with r ms (6 = 0; a = 0.9981) = 1.3 and r ms (6 = —3; a = 
0.9981) = 6.3) and outer radius r ou t = 10, encoded into colors is plotted for representative values 
of tidal charge parameter b = —3.0, 0.0 and inclination of observer #o = 30°, 80°. In the left 
column direct images are ploted, the indirect images are ploted in the central column and the 
composition of direct and indirect images is plotted in the right column. The first two rows of 
images are plotted for the observer inclination #o =30° and the second two rows of images are 
plotted for the observer inclination #o = 80°. Top row images are plotted for b = 0.0, the second 
row images are plotted for b = —3.0, the third row images are plotted for b = 0.0 and bottom row 
images are plotted for b = —3.0. 
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7. Time delay 

For optical effects in vicinity of a black hole, the time delay in case of systems varying 
with time and observed along two different directions due to the light deflection in 
strong gravity can be important. The cordinate time that elapses from the instant 
of photon emission, t e , to the instant of its reception, t , is integrated from the 
Carter equations and reads 



In order to succesfully integrate this formula, one must map all the turning points 
in ji and u motion to correctly set up the signs u sgn and n S gn- 

Suppose that the two light beams, direct and indirect, are emitted at the same 
coordinate time t e . They generally reach the observer at different coordinate times 
^ ir (t 1 G ndir resp.). By time delay we define here the difference At = t 1 G ndir — t dir . 



Fig. 22. The illustration of the impact of tidal charge parameter on the time delay At in case of 
direct and indirect photons emitted from emitter E at coordinate time t e and azimuthal position 
ip e = 7r. They are received at observer O at coordinate times t^ ir (^o ndir resp.). The emittor is on 
circular geodesic in equatorial plane of braneworld Kerr black hole at radial coordinate r = r e . 
The observer is far from the center of the black hole at r = r - Its inclination is = o - 

To demonstrate the impact of the tidal charge b on the time delay we consider the 
following situation (see Figure [22]). Let the isotropicaly radiating monochromatic 
source orbits in the equatorial plane of the braneworld Kerr black hole at radial 
distance r e . It can be switched on and off. When it reaches the azimuthal coordinate 
tp = 7r it is switched on and we compare the coordinate times t dir and t 1 G ndir of 
reception of the photons from the direct and indirect images of the source. 




(107) 



E 



E 
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Fig. 23. The difference ("Time Delay"), At = tjj 



-t^ 11- , between coordinate times of reception 



of direct and indirect geodesies of photons emmited at the same coordinate time t e from the 
azimuthal coordinate cp = 7r is plotted as a function of tidal charge b. Left figure: the inclination 
of the observer is 0q = 20°. Right figure: the inclination of the observer is 6q = 80°. 



The results are demonstrated in the Figure [23j We can directly see that time 
delay At between times of reception of the direct and indirect photons emitted at 
the same instant from the azimuthal position ip = tt increases as the value of the 
tidal charge parameter b goes to higher negative values. When b is fixed, the time 
delay At increases as the value of the inclination decreases. The same effects appear 
for other positions of the radiating spot ((p ^ tt). We can see that the time delay 
At depends strongly on the viewing angle 6q. Therefore, it is extremely important 
to have a system with precisely determined viewing angle. 



8. Optical phenomena related to Sgr A* 

There is an enormously growing evidence that the center of our Galaxy harbors a 
supermassive black hole whose position could be almost surely identified with the 
extremely compact radio source Sgr A*. The chain of arguments seems to be very 
convincing; stars orbiting an unseen mass concentration on elliptical orbits with 
a common focal position, the unseen mass centered on Sgr A* that seems to be 
motionless at the dynamical center of the Galaxy, extremely compact emission of 
the center SSI. Recent measurements of Ghez and collaborators ^01 from the W.M. 
Keck 10 - meter telescopes of a fully unconstrained Keplerian orbit of the short 
period star SO-2 provide the distance Ro = 8.0 ± 0.6 kpc and black hole mass 
M — (4.1 ± 0.6) x 1O 6 M . If the black hole is assumed to be at rest with respect 
to the Milky Way Galaxy (i.e., has no massive companion to induce its motion) 
as argued by Reid the fit can be further constrained to Ro = 8.4 ± 0.4kpc and 
M = (4.5 ± 0.4) x 1O 6 M 

Such a close and huge supermassive black hole could be clearly a very convenient 
object, probably the best one, for testing a wide variety of optical phenomena in 
strong gravity in its vicinity. The time delay of accidents happening behind the 
black hole and observed along two directions could be in principle easily measured. 
We could even expect possibility of black hole silhuette measurements. In this way 
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the influence of the tidal charge could be properly tested and its value estimated, 
because for the Galaxy supermassive black hole we can determine relatively precisely 
the inclination angle of the observer (Solar system), although it is of course very 
close to O ~ 90°. 

For non-rotating , Schwarzchild black holes, the silhuette diameter is given by 
the impact parameter of the photon circular orbit 



D = 2X ph = 6V3M. (108) 

Using the Sgr A* mass estimate M ~ 4.5 x 10 6 M©, we find D c± 55/iarcsec while 
interferometer finges were reported at wavelength of 1.3 mm and fringe spacing 
of 0.00005, comparable with the expected value of D. Shorter wavelengths should 
enable detailed measurements of the black hole silhuette and relatively precise esti- 
mates of the black hole parameters due to very precise knowledge of the inclination 
angle. The angle can be given by the measurement of the Solar system position 
relative to Galaxy plane z e ~ Then 0q ~ 89.9° or more precisely, 0q lies 

between the values of 89.8772° (z Q = ISpc) and 89.9318° (z Q = lOpc). Of course, 
considering the silhuette shape, it is quite enough to take 0q = 90°. 




Fig. 24. Diameter D as a function of braneworld parameter b is plotted for Schwarzchild black 
hole of mass M = 4.5 x 1O 6 M0. Observer is at ro = 8.4kpc lying in the equatorial plane. 



In the case of spherically symmetric black holes, the influence of the tidal charge 
parameter b on the silhuette diameter can be given by the simple formula for impact 
parameter of photon circular orbits that reads 

X ph (b) = 7^= M, (109) 

where 



(110) 
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The resulting dependence of the diameter D(b) is illustrated in Figure [24] • The 
diameter grows slowly with the descending of b; notice that its magnitude is twice 
the pure Schwarzchild value for b = —12.8428. Of course, for rotating black holes the 
silhuette is maximally deformed due to the influence of rotation since the viewing 
angle 0q ~ 90° and is given by calculations and results presented above. Testing of 
the combined spin and tidal charge influence would be possible with measurement 
precision enlarged for 1-order relative to the recently expected state mentioned 
above. Clearly, we can expect that the observational accuracy in near future will be 
high enough to measure the Sgr A* black hole silhuette implying relevant estimates 
of the black hole parameters. 

0.9 

0.8 

0.7 

§0.6 
<l 

0.5 
0.4 
0.3 

-10 -8 -6 -4 -2 

b 

Fig. 25. Comparizon of time delay effect as a function of braneworld parameter b between two 
rotating black holes with rotational parametes a = 0.5 and a = 0.998. For each b the emitter is 
radiating from marginally stable orbit. The relevant values of radii r ms of marginally stable orbits 
are arranged in the Table [6] 



Table 6. Table of relevant values of r ms used in plots oin Fig [25] 



b 




0.0 


-0.5 


-1.0 


-2.0 


-3.0 


-10.0 


r ms (a - 


= 0.5) 


4.24M 


5.05M 


5.73M 


6.88M 


7.85M 


12.88M 


r ms (a = 


0.998) 


1.24M 


3.03M 


3.91M 


5.22M 


6.28M 


11.44M 



Considering the time delay effects, the exact value of 0o is crucial since it plays 
a fundamental role in determining the time delay effect whose scale is given by the 
value of t ~ Isec. We illustrate the influence of the tidal charge on the time delay 
effects at the astrophysically important radii corresponding to marginally stable 
circular geodesies, i.e. in the strong gravity regime, for two representative fixed 
values of black hole spin (see Figure [25] and Table [6]). We can expect importance of 
the regions close to r ms for relevant optical effects due to the idea of the low angular 
momentum accretion in Sgr A* advocated by B. Czerny Clearly, we can see in 
Figure [25] that the time delay effects could be well measurable and the tidal charge 
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influence could be well tested, if the black hole spin is properly estimated. 
9. Conclusions 

One of the most promising ways of estimating influence of hypothetical hidden 
external dimensions, considered in the framework of the braneworld model with 

lOI I 

infinite external dimension as developed by L2JJ , seems to be investigation of the 
optical phenomena caused by the black hole backgrounds. It is so because black 
holes represent the only case when the non-local influence of the bulk space on 
the braneworld spacetime structure can be fully described by a single, braneworld 
parameter called tidal charge, the sign of which can be both positive and negative, 
with the second possibility beeing more realistic one 1 1 \ 

Here, we focused our attention to developing a theoretical background for treat- 
ing the optical phenomena in vicinity of braneworld rotating black holes and bring- 
ing general tendencies of the tidal charge effect in some basical optical phenomena. 

We have shown qualitatively how the braneworld tidal charge affects the basical 
optical phenomena, especially the black-hole silhuette, the accretion disc image with 
the frequency shift of area of the disc radiating at a specific frequency, and the time 
delay between the direct and indirect images of the hot spot orbiting the black hole. 
We have shown that these phenomena could be measured and used to put limits on 
the tidal charge in case of Galaxy Center Sgr A* supermassive black hole. 

We generalized the approaches based on the transfer-func tion method as intro- 
duced and developed in Schwarzchild and Kerr backgrounds H0I 27 | 5 | 42 | 25 | 14 | 17 | 32 ] 
where eq uations of photon motion are solved in terms of the elliptic integrals (see 
| 32 | 23 | 24 ^ p Qr p Ur poses of the present work, the transfer- function method seems to 
be most efficient. Nevertheless, we prepared the ray-tracing method too, since that 
could be usefull in treating other optical phenomena. 

Generally, rising negative value of the tidal charge strenghtens the black hole 
field and suppresses the rotational phenomena, when the black-hole rotation pa- 
rameter is fixed. The magnitude of the optical phenomena grows with decreasing 
of the negatively- valued tidal charge, but the rotation induced asymmetry of the 
phenomena like the black-hole silhuette, or the accretion disc image, decreases. The 
black-hole silhuette is characterized by two parameters, namely the shift of the 
center and ellipticity, that could be in principle measurable in the Galactic Center 
black-hole system Sgr A*, after expected development of observational techniques 
that at present enable measurement of the black hole diameter, not details of the 
shape. The Galaxy center (Sgr A*) seems to be also a promising candidate for test- 
ing the time delay effects both for phenomena related to the accretion disc and 
flares observed there , and for some expected lensing phenomena connected to the 
observed stars orbiting the Sgr ^4* central black hole. 

We have found that observable phenomena could be expected for the time-delay 
effects. Of special interest is comparison of time delays generated for sources in 
vicinity of the Sgr A* black hole (both stars and disc hot spots) and those related 
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to weak lensing of some distant sources I50B3 

Similarly, keeping rotational parameter fixed, the negative tidal charge has ten- 
dency to make the isoradial curve images (both direct and indirect) larger and less 
deformed while the positive tidal charge influence is of opposite character. On the 
other hand, for fixed rotational parameter of the black hole and disc radiating at 
the innermost part above the innermost stable orbit at r = r ms , the negative tidal 
charge restricts the radiating ring image simply because the radius r ms grows with 
decreasing value of braneworld parameter b. Suppresion of the relativistic effects 
can be measurable also in the spectral line profiles generated by the inner hot part 
of the disc radiating at special X-ray line l^SI. 

The optical tests have to be confronted with the data obtained from quasiperi- 
odic oscillations observed in some black-hole systems (micro quasars ^!). The or bital 
resonance model gives good estimates of the black-hole parameters I48 | 46 | 47 [ ^ g 
model has been recently generalized to the case of braneworld Kerr black holes ^1 
It is shown that in the case of microquasar GRS 1915+105 and Galactic Center 
Sgr A* black holes with the negative braneworld parameter b are allowed by the 
observational data^l Detailed modelling of optical phenomena connected to the 
oscillating discs or orbiting (oscillating) hot spots and related resonant phenomena 
between the oscillation modes could be very promising in putting limits on allowed 
values of the tidal charge of the black hole. We plan to elaborate such modelling in 
future. 
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